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Preface: 


ABSTRACT 


The  theory  of  oscillating  thin  arifoils  in  incompressible 
viscous  flow  is  formulated  (see  Refs.  1,  2,  3  for  earlier 
background  material)  and  applied  to  the  calculation  of  steady 
and  unsteady  loads  on  the  family  of  symmetric  Joukowski 
airfoils.  The  theory  is  reduced  to  the  form  of  an  integral 
equation  with  kernel  function  whose  solution  is  obtained  with 
a  modal  expansion  technique  familiar  from  flat  plate  thin 
airfoil  theory.  The  effect  of  viscosity  is  to  change  the  order 
of  the  singularity  in  the  kernel  function  such  that  a  unique 
solution  is  obtained  for  any  cross  sectional  geometry  without 
using  an  auxilliary  uniqueness  criteria  like  the  Kutta 
condition  or  principle  of  minimum  singularity.  The  principal 
unsteady  effect  related  to  thickness  is  to  introduce  an 
explicit  phase  factor  in  the  kernel  function  that  is  propor¬ 
tional  to  the  time  for  vortex  disturbances  to  be  convected 
from  the  load  point  to  the  upwash  point. 

Viscous  thin  airfoil  calculations  for  an  airfoil  with 
sharp  trailing  edge  are  in  agreement  with  the  results  of 
potential  theory  with  Kutta  condition.  The  overall  steady 
loads  calculated  with  either  theory  are  in  general  much  larger 
than  measured  values  and  have  the  wrong  trend  with  increasing 
thickness  ratio.  The  competing  effects  of  thickness  (both 
geometric  and  boundary  layer  displacement)  and  viscosity  for 
an  actual  airfoil  are  such  that  the  measured  loads  are  in 
better  agreement  with  flat  plate  theory  than  with  the 
theoretical  results  for  finite  thickness  airfoils  with  sharp 
trailing  edges. 

Viscous  thin  airfoil  steady  and  unsteady  calculations  for 
an  airfoil  with  elliptic  cross  section  are  in  much  better 
agreement  with  experimental  results.  The  slope  of  the  steady 
state  lift  curve  decreases  with  increasing  thickness  ratio  in 
agreement  with  the  experimental  trend.  The  out  of  phase  loads 
calculated  for  an  oscillating  elliptic,  section  are  significantly 
different  from  either  flat  plate  results  or  those  for  the 
Joukowski  airfoil  with  sharp  trailing  edge.  The  magnitude  of 
the  difference  is  in  quantitative  agreement  with  the  differences 
between  calculations  based  on  Poisso's  integral  equation  and 
recent  experimental  data. 

It  is  concluded  that  viscous  thin  airfoil  theory  is  a 
practical  tool  for  introducing  simultaneously  the  effects  of 
viscosity  and  geometric  thickness  in  two-dimensional  unsteady 
aerodynamic  theory. 
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NOMENCLATURE 


A  coefficients  in  Chebyshev  expansion  of  the  load 

11  distribution,  see  Eq.  (3.7) 

a  reference  length,  radius  of  circle  corresponding 

to  the  Joukowski  airfoils  (quarter  chord  for  flat 
plate  airfoil) 

a,  b  see  Eq.  (3. 18) 


slope  of  the  lift  coefficient  versus  a  curve  at 
zero  lift 
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slope  of  the  leading  edge  moment  coefficient  versus 
a  curve  at  zero  moment 

see  Eqs.  (3.9)  and  (3.12) 

denotes  modified  airfoil  surface,  see  Fig.  2.1 

surface  pressure  coefficient  on  non-lifting  airfoil, 
see  Eq.  (2.2) 

airfoil  chord 

v“  s  +  “ 

exponential  function 

complex  mapping  of  circle  to  Joukowski  airfoil  in 
the  z-plane  and  flat  plate  in  the  -plane,  see 
Eqs.  (3.13)  and  (3.14) 

perturbation  displacement  normal  to  airfoil  surface 
see  Eq.  (3.10) 

enthalpy  field,  non-lifting  airfoil 
perturbation  enthalpy  field 
airfoil  thickness  distribution 

kernel  function,  see  Eqs.  (2.29),  (2.34),  and  (2.47) 
modified  Bessel  function 
integrated  kernel,  see  Eq.  (3.3) 
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WC/2V0,  ,  reduced  frequency 

2wa/v00  ,  modified  reduced  frequency  (reduces  to  k 
for  flat  plate) 

natural  logarithm 

Mach  number 

unit  vector  normal  to  C0  ,  see  Fig.  2.1 

b'/qS 

pressure  field  for  non-lifting  airfoil 

q  J/v» 

surface  speed  on  non-lifting  airfoil 

— 03  t  Reynolds  number  based  on  full  chord 
v 

distance  along  airfoil  surface 
see  Eq .  (2.30)  . 

Chebyshev  polynomial  of  the  first  kind,  Ref.  25, 
p.  774 

time 

unit  vector  tangent  to  C0  ,  see  Fig.  2.1 

Chebyshev  polynomial  of  the  second  kind,  Ref.  25, 
p.  774 

velocity  components  of  non-lifting  flow,  see 
Eq.  (2.19) 

velocity  field  for  non-lifting  airfoil 
perturbation  velocity  field 
free  stream  velocity 
,  surface  upwash 
b0/2av00  (Section  III) 

Cartesian  coordinates  of  leading  and  trailing  edge 
of  flat  plate  airfoil 
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x,  y  Cartesian  coordinates,  see  Fig.  2.1 

AP  ,  see  Eq.  (2 . 23) 
see  Eq.  (2,26) 

W/q  =  |£  +  —  f 

0  q0 

a  angle  of  attack 

Y  0.57721,  Euler's  constant 

AP  see  Eq.  (2.22) 

AcJ>  see  Eq.  (2.22) 

e  eccentricity  parameter  in  Joukowski  mapping,  see 

Eq.  (3.14) 

0(x,y)  see  Eq.  (3.21)  and  Fig.  3.2 

\i  dynamic  viscosity,  p„v 

v  kinematic  viscosity 

v*  see  Eq.  (2.35) 

Pa,  fluid  density 

E  see  Eq.  (3 . 18) 

a  vm  c/4 w  ,  Reynolds  number  referred  to  quarter  chord 

a  v„  a/v  (reduces  to  a  for  flat  plate) 

t  true  thickness  to  chord  ratio 

t*  modified  thickness  ratio 

Tij  viscous  stress  of  non-lifting  flow 

<j>0  potential  of  non-lifting  flow 

4>*  perturbation  potential 

4>a,  <}>r  value  of  <j>0  at  leading  and  trailing  edge  of 
nonlifting  airfoil 
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^  stream  function  of  non-lifting  flow 

q'  perturbation  vorticity 

(jj  radian  frequency  of  simple  harmonic  motion 
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absolute  value  of  quantity  q 
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Laplace  operator 

I.  INTRODUCTION 


The  linearized  inviscid  potential  theory  of  thin  airfoils 
has  served  the  unsteady  aerodynamicist  well  for  many  years. 

In  fact  the  rapid  development  of  modem  high  speed  aircraft 
following  the  second  World  War  would  have  been  severely 
impeded  without  this  powerful  theoretical  tool.  Today,  each 
aircraft  manufacturer  has  its  own  computational  package  for 
carrying  out  aeroelastic  calculations  that  incorporates  a  vast 
array  of  accumulated  theoretical  and  practical  knowledge  on 
the  unsteady  aerodynamic  problem.  When  these  computational 
methods  are  used  by  the  experienced  aerodynamicist  and  with  a 
certain  amount  of  empirical  ’’tailoring”  of  the  unsteady  loads, 
reasonably  good  agreement  with  experimental  measurements  can 
be  obtained.  Unfortunately,  the  most  commonly  used  practical 
index  for  comparison  with  experiment  on  actual  aircraft  has 
been  the  flutter  or  divergence  speed  of  the  wing  -  the  end 
result  of  a  rather  elaborate  aerodynamic  and  aeroelastic 
calculation.  When  theory  and  experiment  do  not  agree  there 
are  many  underlying  assumptions,  both  physical  and  mathematical, 
that  can  be  blamed,  most  of  which  focus  on  some  aspect  of  the 
unsteady  aerodynamic  problem. 

An  obvious  missing  .element  in  thin  airfoil  theory  is  the 
geometric  shape  of  the  airfoil  section.  If  the  flow  is 
incompressible,  inviscid  and  two-dimensional,  the  effect  of 
geometric  shape  is  in  principle  no  problem.  One  can  use  the 
powerful  techniques  of  conformal  mapping  to  obtain  an  exact 
solution  of  the  potential  flow  problem  for  an  airfoil  of  any 
cross-sectional  geometry  that  undergoes  arbitrary  unsteady 
motion.  Unfortunately,  there  is  a  one  parameter  family  of  such 
solutions  for  unsteady  or  steady  flow,  all  of  which  satisfy 
the  boundary  conditions  of  the  potential  flow  problem.  The 
uniqueness  problem  is  even  more  acute  for  the  three-dimensional 
problem,  since  the  spanwise  ’’distribution”  of  the  circulation 
can  be  specified  arbitrarily.  For  the  two-  or  three-dimensional 
potential  flow  problem  the  important  point  is  that  some 
auxilliary  condition  must  be  specified  to  obtain  a  unique 
solution.  In  the  usual  application  of  potential  flow  theory 
it  is  customary  to  apply  the  Kutta  condition  at  the  trailing 
edge  of  the  airfoil. 

The  validity  of  the  Kutta  condition  in  unsteady  aerodynamics 
application  has  been  seriously  questioned  and  debated  in  the 
literature  (e.g.,  see  Refs.  10,  11,  12).  It  is  generally 
recognized  that  the  conventional  thin  airfoil  theory  is  at  best 
unreliable  when  applied  to  hydrofoil  flutter  problems  where  the 
reduced  frequency  is  of  order  unity  or  greater  (Refs.  13,  14, 

15,  16).  The  "smooth”  flow  at  the  trailing  edge  that  is  observed 
in  steady  flow  is  expected  to  break  down  at  some  reduced 
frequency.  In  recent  experimental  work  (Ref.  10)  measurements 
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have  been  made  of  the  surface  pressure  distribution  very  close 
to  the  trailing  edge.  It  was  found  that  locally  large 
departures  from  potential  flow  behavior  can  occur  even  for 
moderately  low  reduced  frequencies. 

The  local  breakdown  of  the  potential  flow  near  the 
trailing  edge  of  an  oscillating  or  fixed  airfoil  should  not  be 
surprising.  If  potential  theory  is  interpreted  as  a  model  of 
the  real  flow  including  the  boundary  layer  and  wake  then  we 
can  safely  argue  that  potential  theory  is  at  its  very  "worst” 
in  the  vicinity  of  the  trailing  edge.  Yet  that  is  precisely 
where  all  of  the  emphasis  is  placed  in  the  usual  application 
of  a  uniqueness  criterion  like  the  Kutta  condition. 
Experimentally,  there  is  no  question  that  the  surface  pressure 
near  the  trailing  edge  of  a  lifting  airfoil  in  steady  flow  is 
not  only  finite  but  almost  fully  recovered  to  the  non- lifting 
value.  However,  if  the  "absence  of  a  trailing  edge  pressure 
singularity"  is  used  as  a  uniqueness  criterion  with  potential 
theory  the  lift  distribution  over  the  entire  airfoil  will  be 
overpredicted  and  the  lift  curve  can  easily  be  2570  too  large. 
Also  the  theoretical  lift  curve  slope  increases,  in  proportion 
to  the  thickness  ratio  while  the  experimental  values  are 
usually  constant  or  decrease  with  increasing  thickness  (see 
Fig.  3.5,  p.  39). 

It  was  shown  by  Pinkerton  (Ref.  17)  that  if  the  measured 
lift  is  used  to  determine  the  unknown  circulation  in  potential 
theory,  then  the  theoretical  and  experimental  pressure 
distributions  are  in  reasonably  good  agreement  (see  Fig.  1.1) 
over  the  forward  957o  of  the  airfoil.  Furthermore,  when  the 
airfoil  shape  was  modified  slightly  (including  a  small  trailing 
edge  bluntness)  to  account  for  boundary  layer  displacement 
effects,  Pinkerton  obtained  near  perfect  agreement  between 
theory  and  experiment  even  to  the  extent  of  eliminating  the 
trailing  edge  singularity.  At  this  point  it  should  be 
emphasized  that  the  only  reason  potential  theory  has  a 
singularity  in  the  first  place  is  because  the  trailing  edge 
is  mathematically  sharp.  It  is  only  for  a  sharp  trailing  edge 
that  the  Kutta  condition  can  be  stated  in  a  mathematically 
precise  way  and  in  general  it  leads  to  serious  errors  even  for 
that  special  geometry. 

The  important  conclusion  from  Pinkerton's  work  is  that  for 
"well  designed"  high  Reynolds  number  airfoils  (i.e.,  no 
separation  over  the  forward  part  of  the  airfoil) ,  inviscid 
potential  flow  theory,  including  geometric  and  boundary  layer 
displacement  effect  is  completely  adequate  to  calculate  the 
detailed  pressure  distribution.  The  only  stipulation  is  that 
one  must  give  up  the  Kutta  condition  or  any  other  uniqueness 
criteria  that  overemphasizes  the  detailed  behavior  of  the  local 
trailing  edge  flow.  The  "fundamental  question"  of  potential 
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flow  theory  still  regains ;  i.e.,  how  can  we  calculate  the 
unknown  circulation?  For  the  unsteady  aerodynamicist,  the 
question  is  particularly  important  since  experimental  measure¬ 
ments  of  the  unsteady  sectional  circulation  are  seldom 
available. 

The  answer  to  the  "fundamental  question’*  of  potential 
theory  must  ultimately  be  found  by  consideration  of  the  direct 
effect  of  viscosity,  the  most  important  missing  physical 
element  of  conventional  thin  airfoil  theory.  It  follows  from 
the  d* Alembert  paradox  (see  Ref.  18,  p.  322  and  405)  that  the 
origin  of  any  net  force  on  a  finite  body  in  an  unbounded  flow 
must  be  cue  to  viscous  effects  no  matter  how  large  the  Reynolds 
number.  The  expression  "inviscid  potential  theory  of  lift"  is 
a  misnomer.  One  must  evaluate  the  lift  in  the  limit  as  the 
?4eynolds  number  tends  to  very  large  values,  a  point  that  was 
emphasized  in  a  recent  paper  of  Sears  (Ref.  19).  The  Kutta 
condition  itself  is  a  statement  about  the  direct  effect  of 
viscosity  although  the  student  of  aerodynamics  rarely  gets 
more  than  a  brief  heuristic  physical  arguement  to  justify  its 
origin. 

It  was  shown  theoretically  by  Shen  and  Crirni  (Ref.  20)  and 
by  direct  calculation  in  the  present  report  that  the  Kutta 
condition  is  the  correct  viscous  uniqueness  criterion  in  the 
high  Reynolds  number  limit  for  airfoils  with  a  mathematically 
sharp  trailing  edge.  The  proof  of  uniqueness  of  the  thin 
airfoil  problem  using  Oseen's  theory  as  in  Ref.  20,  or  viscous 
thin  airfoil  theory  as  in  this  report,  should  be  included  in 
any  modern  text  bock  or  college  level  course  on  the  theory  of 
lift.  The  flat  plate  the  ry  (see  Section  II. C)  is  straight¬ 
forward  and  permits  cne  to  see  the  explicit  connection  between 
the  effect  of  viscosity  in  the  Stokes  layer  and  the  lift  on  an 
airfoil. 

The  practical  limitation  of  the  flat  plate  theory  is  that 
direct  viscous  effects  are  "grossly"  overexaggerated  because  of 
the  mathematically  sharp  trailing  edge.  Recall  from  Pinkerton's 
work  (Ref.  17)  that  the  seme  criticism  can  be  leveled  at  the 
Kutta  condition  (also  see  Section  III.C).  The  trailing  edge 
Reynolds  number  (i.e.,  based  on  the  edge  thickness)  of  typical 
airfoils  Is  of  the  order  of  IQ3  or  10' .  If  the  Reynolds  number 
is  based  cn  the  displacement  thickness  at  the  trailing  edge  the 
numbers  are  more  like  IQ*  or  10s.  The  geometric  edge  is 
totally  immersed  in  a  turbulent  shear  flow  and  the  direct 
effect  of  viscosity  near  the  trailing  edge  is  diminished.  The 
details  of  the  trailing  edge  flow  practically  defy  analysis  and 
sericus  attempts  to  analyze  the  problem  have  been  restricted 
to  airfoils  with  sharp  trailing  edges  (see  Refs.  21,  22,  23). 

It  is  a  basic  premise  of  the  present  author  that  the  ultra- 
detail  of  the  trailing  edge  flow  cannot  be  of  serious  practical 
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importance  to  the  unsteady  aerodynamicist .  The  arguments  to 
support  this  claim  are  based  on  a  mixture  of  experimental  and 
theoretical  evidence.  First  of  all  the  overall  lift  and  moment 
(and  higher  order  moments  required  for  aeroelastic  analysis) 
of  an  airfoil  at  high  Reynolds  number  ( > 106)  are  "linear" 
functions  of  the  angle  of  attack.  The  trailing  edge  flow  is 
very  nonlinear  (separated  and  turbulent)  and  yet  the  non¬ 
linearity  does  not  appear  in  the  overall  loads  until  the  stall 
speed  is  approached.  On  the  other  hand,  a  nonlinear  leading 
edge  flow  has  a  profound  nonlinear  effect  on  the  overall  loads 
(see  Section  III.C,  Fig.  3.7). 

Other  evidence  to  support  our  claim  of  "trailing  edge 
independence"  may  be.  inferred  from  the  recent  unsteady  airfoil 
experiments  of  Davis  and  Malcolm  (Ref.  24) .  All  of  their 
unsteady  surface  pressure  measurements  exhibit  a  linear  simple 
harmonic  relation  to  the  forced  motion  of  the  airfoil.  No 
second  harmonic  content  of  any  magnitude  was  detected  in  the 
measurements.  Also,  the  measured  unsteady  loads  are  virtually 
zero  over  the  aft  10%  of  the  chord  for  reduced  frequencies  up 
to  0.3.  Therefore,  the  trailing  edge  loads,  whether  they  are 
linear  or  nonlinear,  do  not  contribute  to  the  overall  load 
coefficients  required  in  aeroelastic  analyses. 

The  effect  of  viscosity  at  the  trailing  edge  is  de- 
emphasized  in  the  present  work  by  considering  an  effective 
trailing  edge  bluntness,  thus  eliminating  the  associated 
potential  flow  singularity  as  in  the  work  of  Pinkerton.  The 
direct  effect  of  viscosity  over  the  entire  geometric  cross 
section  corrected  for  boundary  layer  displacement  is  the 
physical  origin  of  the  circulatory  lift.  The  calculation  of 
the  unsteady  loads,  including  geometric  and  viscous  effects, 
is  the  primary  aim  of  viscous  thin  airfoil  theory.  The  basic 
mathematical  concepts  are  formulated  in  detail  in  Section  II. 

An  integral  equation  is  derived  that  relates  the  loads  and 
upwash  via  a  kernel  function  that  includes  the  essential 
effects  of  geometry  and  viscosity.  Numerical  calculations  of 
the  loads  on  fixed  and  oscillating  airfoils  of  the  Joukowski 
family  are  compared  with  experimental  data  in  Section  III. 


II.  VISCOUS  THIN  AIRFOIL  THEORY 

A.  Basic  Concepts  and  Equations 

Consider  a  thin  two-dimensional  symmetric  airfoil  in  an 
incompressible  flow  at  zero  angle  of  attack  (Fig.  2.1a).  For 
very  high  Reynolds  number  the  nonlifting  pressure  distribution 
may  be  calculated  quite  accurately  with  potential  flow  theory 
(e.g.,  see  Fig.  1.1).  If  the  thin  boundary  layer  displacement 
thickness  is  added  to  the  geometric  thickness  (dashed  line  in 
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Figure  2.1.  -  Geometry  of  the  viscous  thin  airfoil  problem  and  removal  of  the 
symmetric  boundary  layer. 


Fig.  2.1a)  even  greater  accuracy  may  be  obtained.  Alter¬ 
natively,  given  the  measured  surface  pressure  distribution 
on  an  actual  body,  one  can  solve  the  inverse  potential  flow 
problem  and  determine  the  precise  shape  of  a  ''modified  body" 
that  will  induce  such  a  pressure  distribution.  Except 
perhaps  in  the  immediate  vicinity  of  the  trailing  edge, 
the  difference  between  the  modified  and  actual  body  shapes 
will  be  of  the  order  of  the  boundary  layer  displacement 
thickness . 

The  nonlifting  potential  flow  so  constructed,  may  be 
interpreted  in  a  slightly  different  way.  The  modified  airfoil 
in  Fig.  2.1a  may  be  thought  of  as  a  real  airfoil  whose  upper 
and  lower  surfaces  move  at  precisely  the  potential  flow 
speed,  q0  ,  calculated  from  the  inverse  Bernoulli  relation 
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(2.1) 


with  the  surface  pressure  coefficient 


Po  ~  P 
ho  v2 

2  Hqo  vco 


b.Q  -  ho, 


(2.2) 


The  potential  flow  solution  for  the  modified  airfoil  may  then 
be  interpreted  as  an  exact  viscous  solution  of  the  Navier- 
Stokes  equations;  i.e., 


div  v0  =  0 


(2.3) 


•  gradv0  +  grad  h0  =  v  V2  vc 


(2.4) 


n  •  v0  =  0 


t  *  v0  =  q0 


on  C0 


(2.5) 


v„  v  V,v, 
'0  00 


at  00 


(2.6) 


The  vorticity  is  zero  everywhere  so  that 
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grad  4>0 


(2.7) 


The  viscous  stress  is  given  by  the  expression 


2y 


32  4>0 

Sx1  3x^ 


(2.8) 


It  is  not  zero  but  its  divergence  is  zero  so  that  v0  ,  h0  are 
also  solutions  of  the  inviscid  Euler  equations. 

The  purpose  of  the  foregoing  construction,  is  to  remove 
the  thin  turbulent  boundary  layer  and  wake  on  the  actual  airfoil. 
The  geometric  airfoil  shape  has  been  modified  slightly  to  account 
for  the  high  Reynolds  number  boundary  layer  displacement 
thickness.  We  are  left  with  an  exact  viscous-potential  flow 
over  a  modified  body  C0  with  a  moving  surface  and  a  surface 
pressure  distribution  that  is  as  close  to  the  actual  measured 
pressure  distribution  as  we  choose  to  calculate.  Finally,  we 
emphasize  the  point  that  the  constructed  equivalent  airfoil  and 
potential  flow  will  represent  the  actual  airfoil  and  real  flow 
’’best"  at  high  Reynolds  number.  It  is  tacitly  assumed,  of 
course,  that  the  actual  airfoil  is  "designed"  for  high  Reynolds 
number  operation  with  an  attached  boundary  layer  flow,  except 
perhaps  near  the  trailing  edge. 

We  now  pose  the  fundamental  question  of  viscous  thin 
airfoil  theory.  What  load  distribution  is  induced  on  the 
"modified"  airfoil  when  it  is  perturbed  from  its  nonlifting 
configuration?  (See  Fig.  2.1b.)  The  surface  displacement 
and/or  velocity  is  assumed  to  be  small  but  may  otherwise  be 
steady  or  unsteady.  The  modified  airfoil  surface  is  assumed  to 
undergo  the  same  deformation  as  the  actual  airfoil  surface; 
i.e.,  the  boundary  layer  thickness  is  assumed  to  be  "frozen". 

The  perturbed  flow  (denoted  by  a  primed  variable)  satisfies  the 
following  set  of  linearized  viscous  equations: 


divv’  =  0 


(2.9) 


Do  v' 
Dt 


+  v'  •  grad  v0  +  sradh' 


v  V2  v' 


(2.10) 


where 
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(2.11) 


v  - 


h'  =  h  -  h0 


and 

=  _L  +  vQ  •  grad  (2.12) 

T\+-  's*.  O  ° 


Also,  v'  must  vanish  at  infinity  and  satisfy  viscous 
boundary  conditions  on  the  modified  airfoil  surface;  i.e., 

v  =  v0  +  v'  =  Vq  on  C  (2.13) 


(see  Fig.  2.1b),  where  v0  is  the  specified  surface  velocity. 
For  the  lifting  problem  we  assume  that  the  perturbation 
displacement  and  velocity  are  everywhere  normal  to  the 
unperturbed  surface,  C0  .  Then  the  projection  of  the  boundary 
conditions  onto  C0  (see  Fig.  2.1b)  yields  the  following  two 
conditions : 


v' 


on  C0 


(2.14) 

(2.15) 


where  f  is  the  perturbation  displacement  normal  to  the 
unperturbed  surface.  The  first  boundary  condition  :s 
recognised  as  the  familiar  inviscid  requirement  of  no  flow 
through  the  surface.  The  second  condition  is  a  "viscous" 
no-slip  condition  on  the  "perturbation  flow".  The  fundamental 
problem  of  viscous  thin  airfoil  theory  is  to  solve  Eqs.  (2.9) 
and  (2.10)  in  the  limit  as  v  0  ,  subject  to  the  boundary 
conditions  (2.14)  and  (2,15)  and  the  requirement  of  a 
vanishing  perturbation  flow  at  infinity. 


B.  The  Associated  Potential  Flow  Problem 

If  the  viscosity  is  set  equal  to  zero  in  Eq.  (2.10)  and 
the  viscous  boundary  condition  Eq.  (2.15)  is  dropped,  we 
recover  the  familiar  "inviscid"  linearized  potential  flow 
problem  for  the  load  distribution.  Introduce  a  perturbation 
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velocity  potential  such  that 


grad  <J> ' 


(2.16) 


then  the  "inviscid"  problem  becomes 


V2<}>'  =  0 


Roll 

Dt 


=  W(x,  t) 


on  Cr 


|v'  I  v  0  at  «> 


(2.17) 


In  addition  to  the  boundary  conditions  on  the  surface  and  at 
infinity,  we  must  impose  one  additional  and  all  important 
condition  to  obtain  a  unique  solution  of  the  potential  flow 
problem.  For  the  moment  we  leave  this  uniqueness  condition 
open. 

The  potential  flow  problem  can  of  course  be  solved 
exactly  by  conformal  mapping  techniques.  Here  we  derive  an 
integral  equation  from  which  we  can  subsequently  formulate  the 
viscous  problem  for  a  body  with  thickness.  In  the  remainder 
of  this  report  we  assume  that  all  perturbation  variables  have 
a  simple  harmonic  time  dependence;  i.e., 


q(x,t)  =  q(x)  elut 


(2.18) 


where  w  is  the  radian  frequency.  Next  we  introduce  a 
conformal  coordinate  transformation  in  Eqs.  2.17  that  maps 
the  airfoil  onto  a  '‘slit'1  on  the  real  axis.  If  $0,  are 
the  conjugate  potential  and  streamf unction  corresponding  to 
the  thickness  induced  flow  field,  i.e., 


(2.19) 
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then  $0,  are  suitable  mapping  functions.  In  <J>0,  tyo 
coordinates  the  potential  flow  problem  becomes  (see  Fig.  2.2) 


Let 


Then 


V2 

3  2  <{> ' 

+ 32 

3*! 

3  'I'o 

3<j> ' 
3*0 

^o-0± 

“  W/q0 

«  yr 

ML  + 

3<j>o 

if*' 

qo 

qo 

-  -  p 

| grad  <j> '  | 

'V;  0 

A(j>  - 

i 

o 

o 

-e- 

-e- 

-  4>'  (4»0 

AP  * 

P(*o>0") 

-  P(4>o.0 

1M 

3*o 

,  iic  .  , 

+  -y  A(j) 

qo 

=  -  AP 

<  4>0  <  ^ 


B 


Ah' 

qo 


and 


i  o 


A4>  =  -  f  Jz?(<l>o)  exp 


-i.  f° 

*o 


qo 


The  solution  of  Eq.  (2.20)  may  be  expressed  in  terms  of 
potential  doublets  on  the  real  axis;  i.e., 


(2.20) 


(2.21) 


(2.22) 


(2.23) 


(2.24) 
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♦  ’  -  -  ft  (2.25) 

♦a 

where 

&  -  [<4>0  -  *i)2  +  *§)%  (2.26) 

Now  for  t|?0  0 

00 

ft  ■  h  /*♦<♦*>  ft  <2-27> 

Substitute  Eq.  (2.24)  into  Eq.  (2.27),  interchange  the  order 
of  integration  and  take  the  limit  as  \J|q  -►  0  .  We  obtain  the 
following  integral  equation  for  the  load  distribution: 


4>b 

i  /  jgtM)  K«o-*oi  *o)<Mo  -  -Wl$o)  (2.28) 

*A 

where 

K(*»i*i)  =  ft  /  ds  “pf-1”  /  if)  ft  tal*o-=l 

0  \  *f»o  / 

(2.29) 

and  the  slash  through  the  integral  indicates  a  Cauchy  principle 
value.  The  inviscid  integral  equation  (2.28)  and  kernel 
function  K(<j>0;<j>o)  will  play  central  roles  in  the  theoretical 
development  of  viscous  thin  airfoil  theory. 

The  non-uniqueness  of  the  potential  flow  problem  is 
readily  apparent  from  the  integral  formulation.  The  kernel 
function  (2.29)  has  a  Cauchy  singularity  at  4>0  =  0  .  It  follows 
that  the  integral  equation  has  an  eigensolution  [see  Eq.  (2.39)]. 
Some  auxilliary  condition  must  be  specified  if  inviscid 


13 


potential  theory  alone  is  to  be  useful.  The  alternative  is  to 
develop  a  viscous  theory  or  to  derive  an  appropriate  auxilliary 
condition  from  viscous  theory.  That  is  the  primary  objective 
of  this  report. 

Even  though  the  solution  of  the  potential  theory  integral 
equation  is  not  unique,  the  essential  effects  of  geometric  and 
boundary  layer  displacement  thickness  are  explicitly  brought 
out.  First  of  all  the  upwash  and  the  load  distribution  are 
modified  by  the  local  surface  speed.  The  unknown  load  function 

is  actually  the  pressure  coefficient  based  on  the 
"local"  dynamic  pressure.  The  most  important  effect  of  thick¬ 
ness  on  the  unsteady  problem  is  to  induce  the  variable  phase 
factor  in  the  kernel  function.  The  quantity 

+  s 

T  =  f  dt/q2  (2.30) 

♦o 


is  physically  the  time  for  a  vortex  disturbance  to  be  convected 
from  the  load  point  on  the  airfoil  surface  to  the  point  where 
the  upwash  is  observed.  For  the  flat  plate  airfoil,  T  is 
proportional  to  the  distance  between  the  load  and  upwash  points 
since  the  local  convection  velocity  is  the  constant  free  stream 
velocity.  For  a  body  with  finite  thickness  the  variable  phase 
becomes  an  important  factor  in  the  calculation  of  the  out  of 
phase  component  of  the  load.  Because  of  this  factor,  we  also 
note  that  the  kernel  (2.29)  is  not  translation  invariant  as 
it  is  for  a  flat  plate.  In  the  integral  equation,  the  unsteady 
kernel  K(4>o  “  $6;  <5>o)  depends  on  the  explicit  location  of  the 
load,  4>o  »  as  well  as  the  difference,  4>o  ~  »  between  the 

load  and  upwash  points . 


C.  Zero  Thickness  Viscous  Theory 

We  return  now  to  the  viscous  perturbation  problem  of 
Section  A.  For  zero  thickness  and  simple  harmonic  motion  the 
perturbation  problem  becomes 

divv’  *  0 


Deo  v’ 

Dx 


+  gradh' 


v  V2  v' 
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■t  , 

J  •  v' 


4*  ■+•  I 

1  •  V 


4. 

on  y  =  0 


XA  <  x  <  XB 


|v‘ |  ^  0 


at  ® 


(2.31) 


where 


Deo  a  .  . 

—  =  Veo  —  +  rw 

Dx  3x 


(2.32) 


is  the  simple  harmonic  convective  derivative.  In  a  previous 
report  (Ref.  2)  an  exact  integral  equation  was  derived  for  the 
zero  thickness  viscous  problem.  Here  we  express  the  kernel  in 
the  notation  of  Section  B.  With 


q0  -  v0 


x  -  $0/Voo  „  y  =  v0/v« 


(2.33) 


the  exact  flat  plate  kernel  function  is 


K(*0)  -  J-  /  d*;  «p(- 


where 


X  3^7  An^o  "  *ol  +  exP 


* 

^$0  ~  $0  ^  ^  ^  1 4>p  -  <j>o  1  ^ 


(2.34) 


v*  =  v  /  1  +  2i 


(2.35) 
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and  K0(z)  is  the  modified  Bessel  function  (Kef.  25).  The 
solution  of  Eq.  (2.2S;  with  the  flat  plate  viscous  kernel, 

Eq.  (2.34),  was  the  subject  of  Ref.  2,  where  we  also  discussed 
the  relationship  of  viscous  thin  airfoil  theory  to  the  theory 
of  lift  based  on  Oseen's  equation  (see  Ref.  20).  The  main 
points  are  summarized  below. 

For  any  value  of  the  viscosity,  however  small,  the 
Cauchy  singularity  in  the  inviscid  kernel  is  reduced  to  a 
symmetric  logarithmic  singularity  as  shown  in  Fig.  2.3.  With 
a  logarithmic  singularity  there  are  no  eigensolutions  of  the 
integral  equation  (Ref.  26).  The  solution  of  the  viscous  flat 
plate  airfoil  equation  is  unique.  As  v  -»■  0  or  Reynolds 
number  Re  tends  to  infinity,  the  solution  tends  to  the 
potential  flow  result  with  the  Kutta  condition  applied  at  the 
trailing  edge.  The  results  shown  in  Fig.  2.4  illustrate  this 
point  for  the  f^at  plate  at  constant  angle  of  attack.  For 
high  Reynolds  number,  the  difference  between  the  viscous  and 
inviscid  solutions  is  of  order  1//Re  except  near  the  trailing 
edge,  a  result  due  to  Shen  and  Crimi  (Ref.  20).  The  calculated 
deviations  from  potential  flow  behavior  near  the  trailing  edge 
are  not  large  enough  to  account  for  measured  local  deviations 
reported  in  Ref.  10.  A  summary  of  these  detailed  comparisons 
may  be  found  in  Ref.  2. 

Far  downstream,  the  effect  of  viscosity  is  to  cause  an 
exponential  decay  of  the  oscillating  wake.  For  ordinary 
aerodynamic  applications  (v^  »  /wv  )  the  wake  decay  is  of 
negligible  importance  and  the  complex  kinematic  viscosity  v* 
in  the  kernel  [see  Eq.  (2.35)]  may  be  replaced  by  v  ,  thus 
simplifying  the  actual  evaluation  of  the  kernel.  The  direct 
effect  of  viscosity  is  quasi-steady. 

Some  additional  features  of  the  flat  plate  viscous  theory 
(not  reported  in  Ref.  2)  can  be  brought  out  by  comparing  the 
integral  formulation  and  solution  for  very  low  Reynolds  number 
(Stokes  flow)  with  the  corresponding  inviscid  results.  For 
simplicity  consider  the  special  case  of  steady  flow  over  a 
flat  plate  airfoil  at  constant  angle  of  attack,  a  .  Introduce 
the  dimensionless  coordinate  x  «  2§0/vm  c  where  c  is  the 
airfoil  chord.  Then  the  low  Reynolds  number  and  infinite 
Reynolds  number  integral  equations  and  solutions  are 
respectively: 

For  o  =  Veo  c/4v  •»  0  (Stokes  Flow) : 

1 

~  /  J^tS)  [  £n|x  -  5 1  +  in  a  +  y  +  1  ]  d£  =  -§ 

-1 
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(2.36) 


Figure  2,3.  -  Kernel  function  of  tho  stondy  viscous  airfoil  equation. 


Figure  2.4.  -  Convergence  of  solution  of  the  viscous  airfoil  equation  for  a  flat  plate  at 
constant  angle  of  attack. 


where  y  s  0,57721  is  Euler’s  constant.  The  "unique"* 
solution  is : 

•£*«  - - s-2 -  (2'37) 

a  to  |  +  y  +  1  -  x 


where  A  is  an  arbitrary  constant,  since  1/A  -  xr  is 
an  eigensolution  of  Eq.  (2.38). 

The  load  distribution  in  the  Stokes  flow  limit  is 
precisely  the  same  form  as  the  "eigensolution"  of  inviscid 
theory.  Although  we  have  used  high  Reynolds  number  arguments 
to  formulate  the  viscou  thin  airfoil  problem  there  is 
nothing  special  about  t..a  formal  application  of  the  theory 
for  low  Reynolds  number.  In  fact  we  expect  a  continuous 
dependence  of  the  solution  on  Reynolds  number.  It  follows 
that  for  a  flat  plate  airfoil  with  mathematically  sharp 
leading  and  trailing  edges  there  must  be  a  square  root 
singularity  at  both  edges ,  even  though  the  flow  is  viscous . 

It  is  only  at  high  Reynolds  number  that  we  see  a  tendency  to 
weaken  the  trailing  edge  singularity.  In  this  sense  the 
mathematical  statement  of  the  Kutta  condition  is  a  limiting 
condition  for  high  Reynolds  number  flow. 

The  most  revealing  feature  of  the  foregoing  comparison 
is  that  the  ultimate  origin  of  lift  for  any  finite  value  of  the 
Reynolds  number  is  through  the  direct  action  of  viscosity  in 


*For  a  general  proof  of  uniqueness  of  Stokes  flow  see 
Batchelor  (R.ef.  18,  p.  227). 


the  Stokes  boundary  layer.  The  eigensolution  of  so-called 
"inviscid"  theory  is  actually  the  last  remnant  of  the  viscous 
effect,  as  it  becomes  more  and  more  localized  near  the 
airfoil  surface  in  the  high  Reynolds  number  limit.  The 
auxilliary  uniqueness  statement  (e.g.,  the  Kutta  condition, 
or  principle  of  minimum  singularity)  that  must  be  imposed  on 
the  family  of  potential  flow  solutions  (see  Section  II. B,  p.  9) 
is  a  statement  about  the  "direct"  effect  of  viscosity.  The 
expression  "inviscid"  potential  theory  of  lift  is  a  misnomer 
since  it  is  precisely  the  "singular"  action  of  viscosity  in 
the  high  Reynolds  number  limit  that  is  ultimately  responsible 
for  the  lift.  This  point  was  emphasized  in  a  recent  paper  by 
Sears  (Ref.  19). 


D.  Viscous  Theory  of  Thin  Airfoils 

We  now  formulate  the  fundamental  integral  equation  of 
viscous  thin  airfoil  theory.  The  main  result  follows  by 
comparison  of  the  "exact"  potential  flow  kernel  that  includes 
the  effect  of  thickness  [Eq.  (2.29)]  with  the  "exact"  viscous 
flow  kernel  for  the  zero  thickness  airfoil  [Eq.  (2.34)].  The 
key  observation  is  that  for  high  Reynolds  number  and  moderate 
reduced  frequency  the  important  viscous  effect  on  the  kernel 
is  to  reduce  the  order  and  modify  the  potential  flow  singularity 
over  a  distance  of  0  (c//Re  )  [see  Eq.  (2.34)  and  Fig.  2.2]. 

For  Voo  »  /unT  ,  a  condition  that  is  always  satisfied  in  ordinary 
aerodynamic  applications,  the  effect  of  viscosity  is  quasi¬ 
steady;  i.e.,  v*  =  \)  in  Eq.  (2.35).  Also  the  viscous  effect 
is  additive  in  the  kernel  [see  Eq.  (2.34)]. 


Now  take  the  curl  of  the  perturbation  momentum  equation 


(2.10) 

and  transform  into  <j>0,  \p0  coordinates. 

The  result  is 

where 

v  V2  O'  -  ( —  +  iy-  =  0 

\a*o  9o  / 

(2.40) 

O'  -  k  •  curl v' 

(2.41) 

is  the 

perturbation  vorticity.  Now  let 

O’  =  exp (<f>0/2v)  Q(<j>0,i|>0) 

(2.42) 

Then 
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u 


K> 


'  J 


V2  Q  -  (i  +  4i^Q  =  0  (2.43) 

and  for  q0  »  *4Sv  ,  the  quasi-steady  approximation,  the 
fundamental  solution  of  (2.43)  is  the  modified  Bessel  function 
so  that 


where 


8'  =  exp(<j>0/2v)  K0(^/2v) 


(2.44) 


&  -  <*§  + 


(2.45) 


The  local  "quasi-steady"  viscous  solution  for  a  body  with 
finite  thickness  is  precisely  of  the  same  fora  in  $0>  \J>0 
coordinates  as  the  corresponding  flat  plate  expression  (see 
Ref.  27).  By  analogy,  the  integral  equation  and  kernel  of 
viscous  thin  airfoil  theory  can  be  expressed  as  follows: 


_L 

2tt 


oSt'J’o)  K(<f>o  "  4> o »  4>o)  d  4>o 


where 


-3T(*0>  (2.46) 


K(<J>qJ  *£) 


00 

■—  f  ds  exp 
9<p0  *7 

o 


x 


s|  +  exp[(<j>0  -  s)/2v]  K0 c 1 4>0  -  s|/2v) 


(2.47) 

The  last  two  equations  are  the  fundamental  theoretical  results 
of  this  report.  Note  that  the  kernel  has  a  symmetric  logarith¬ 
mic  singularity  so  that  all  of  the  integrals  are  to  be 
interpreted  in  the  ordinary  sense.  The  weak  singularity  also 
implies  that  the  solution  of  (2,46)  is  unique.  Numerical 
calculations  based  on  these  equations  are  summarized  in 
Section  III. 
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III.  NUMERICAL  RESULTS 
A.  Method  of  Solution 

To  solve  the  integral  equation  (2.46)  we  first  introduce 
the  dimensionless  variable 


v  -  4*0 

2a  Vo, 


(3.1) 


where  a  is  a  length  to  be  chosen  such  that 

4>o  “  4> a  when  x  =  -1 

(3.2) 

4>0  =  <j)g  when  x  =  1 

then  the  integral  equation  and  kernel  can  be  written  in  the 
dimensionless  form 


X 

i  f  SHj)  K(x-y;  y)  dy  =  -  (3-3) 

-1 


where 


K(x;  y)  -  —  Kj  (x;  y) 


(3.4) 


and 


Kj(x;y)  =  J  d£  exp | 


y  +  t 

-  ik  J  dx/Qc 


3x 


£&n | x  -  £  |  +  exp[3(x- ?)]  K0(5|x- 5|)] 


(3.5) 


Also 
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Qj  «  qo/v~ 

£  =  2w  a/vTO 


a 


Vos  a 
v 


^  = 


1 


ar- 


3f  ,  ik  r 

w  +  q7  £ 


(3.6) 


The  last  expression,  the  local  dimensionless  surface  upwash, 
must  be  expressed  in  terms  of  the  normalized  potential  coordi¬ 
nate,  x .  The  parameter,  s  ,  is  the  distance  along  the  airfoil 
surface.  The  variable  x  should  not  be  confused  with  the 
streamwise  Cartesian  coordinate.  We  refer  to  Kj(x;y)  as  the 
integrated  kernel. 

Next  we  expand  the  load  function  Jzi^x)  in  a  series  of 
Chebyshev  polynomials ;  i.e., 


&U) 


N 


£ 

n  =  0 


T_  (x) 
A  ..-IV 1 


(3.7) 


where  Tn(x)  is  the  Chebyshev  polynomial  of  the  first  kind 
(Ref.  25).  Substitute  (3.7)  into  (3.3),  multiply  by  /I  -  xz  Um(x) 
(Chebyshev  polynomial  of  the  second  kind)  and  integrate  from 
-1  to  1.  The  end  result  is  a  matrix  equation  for  the  unknown 
load  coefficients,  ^  ;  i.e.  , 

N 

£  Cmn  An  =  So  ”  '  0,  (3.8) 

n  =  0 


where 
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_  (m  +  1)  f  Tn(y>  f  Tm+1(x) 


mn 


7T 


/  -ILLt  dy  f  -^±1=4-  Kj  (x  -  y;  y) 
<  /I  -  yz  <  /I  -  xz 


/  l  -  yT 

-1  y  -1 


dx 


and 


(3.9) 


Jm 

I  /  um(x)  STix)  dx  (3.10) 


-1 


If  we  compare  the  last,  three  equations  with  the  corresponding 
flat  plate  results  [see  Eqs.  (5.2),  (5.3)  and  (5.9)  of  Ref.  2] 
the  main  difference  is  that  the  unsteady  integrated  kernel  in 
(3.9)  is  an  explicit  function  of  the  load  point  as  well  as  the 
difference  between  the  load  and  upwash  points  -  a  direct 
consequence  of  the  variable  surface  speed.  It  is  a  curious 
fact  that  the  steady  state  kernel  (k  =  0)  ,  including  thickness, 
is  independent  of  the  explicit  load  position.  The  spectral 
coefficients  (3.9)  of  the  integrated  kernel  are  not  only 
functions  of  the  reduced  frequency  and  Reynolds  number,  as  in 
the  flat  plate  case,  but  they  are  also  functions  of  the  cross- 
sectional  geometry.  The  coefficients  depend  on  the  geometry 
both  through  the  mapping  functions  (<|>0,  i}»0)  and  through  the 
unsteady  convective  phase  factor  in  the  kernel.  The  actual 
numerical  evaluation  of  the  spectral  coefficient  is  simplified 
with  the  trigonometric  substitution 

x  =  cos  0  ,  y  =  cos  4>  (3.11) 


in  which  case 

TT  IT 

^mn  ~  J  cos  n  d-4>  J  cos[  (m  +  1)  8]  d6 

o  o 


x  Kj  (cos  8  -  cos  4> ;  cos  <J>)  (3.12) 


The  kernel  is  first  evaluated  for  a  discrete  set  of  load  points 
as  a  function  of  the  difference  argument.  Then  the  's 
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are  evaluated  by  successive  quadratures  in  the  order  indicated 
in  Eq.  (3.12).  Since  the  integrated  kernel  Kx  is  non¬ 
singular,  the  numerical  integration  poses  no  formal  difficulty. 
The  time  consuming  part  of  the  operation  is  the  detailed 
evaluation  of  the  kernel. 


B.  The  Kernel  Function 

It  is  a  remarkable  feature  of  two-dimensional  viscous 
incompressible  thin-airfoil  theory  that  a  kernel  function  for 
an  oscillating  airfoil  with  finite  thickness  can  be  derived, 
and  furthermore,  that  in  the  mean  flow  coordinate  system  it 
has  the  same  structure  as  the  corresponding  flat  plate  result. 
The  steady  flow  induced  by  the  geometric  and  boundary  layer 
displacement  thickness  plays  a  most  important  role  in  the 
structure  of  the  kernel. 


In  the  remainder  of  this  report  we  consider  the  thickness 
induced  flow  for  the  symmetric  Joukowski  family  of  airfoils. 
The  complex  thickness  potential  and  mapping  are  given  below. 


(3.13) 
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On  the  surface  of  the  airfoil 


w 


a  e 


ie 


(3.15) 


<!>0  =  2vxa  cos  9  =  <f>Q-  +  ii|>0  (3.16) 


and 


e  + 


(3.17) 


From  (3.16)  we  see  that  the  proper  length  scale  in  the 
dimensionless  potential  variable,  x,  in  (3.1)  is  the  radius 
of  the  airfoil  circle  in  the  w-plane.  The  parameters  a,  b 
are  related  to  the  airfoil  chord,  c  ,  thickness  ratio,  x*  , 
and  eccentricity,  s  ,  by  the  following  formulae: 
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C  1  +  T* 

4  E 


(3.18) 


For  the  elliptic  cross  section  (e  =  0) ,  t*  is  the  true 
thickness  ratio.  For  e  ^  0,  t*  is  slightly  less  than  the 
true  thickness  ratio  and  is  calculated  with  the  thickness 
corresponding  to  the  point  w  =  ia  in  the  circle  plane.  The 
square  of  the  speed  on  the  stagnation  stream  line  is  given 
by  [see  Eq.  (3.6)] : 


q2  =  [1  “  a2/w2|^ 

°  |F* (w)|* 

4  0 

w  =  a  e  0  <  0  <  it  on  the  body 

w  =  at  1  <  t  <  °°  in  the  wake  (3.19) 

In  the  kernel,  Eq.  (3.5),  we  have, 

x  =  cos  9  ,  0  <  6  <  TT 

x  ”  hit  +  l/t)  ,  1  <  t  <  «•  (3.20) 


and  the  dimensionless  convective  phase  factor 


6<x’y)  =  /  | 
x 


-1  <  X  <  1 

x  <  y  <  00 


(3.21) 
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Also  note  that  the  reduced  frequency  k  and  Reynolds  number 
'  a  are  based  on  the  length  scale  a  in  the  kernel.  For  the 
flat  plate  a  =  c/4  ,  the  reduced  frequency  scales  with  the 
semi- chord  and  the  Reynolds  number,  Re  ,  is  referred  to  the 
full  chord;  i.e.  , 


2Voo 

Re  =  (3.22) 


Typical  plots  of  Qq  are  given  in  Fig.  3.1a  and  b  for 
the  elliptic  cross  section  (t  =  12%)  and  the  Joukowski  airfoil 
with  cusped  trailing  edge  (t  =  157«) .  The  variable  surface 
speed  is  integrated  numerically  to  yield  the  typical  convective 
phase  factor  illustrated  in  Fig.  3.2  [see  Eq.  (3.21)].  The 
results  shown  are  for  an  elliptic  cross  section  approximately 
10%  thick.  The  load  point  is  at  the  leading  edge  and  the 
vertical  difference  between  the  dashed  and  solid  lines  is 
proportional  to  the  difference  in  convection  times  for  a 
vortex  disturbance  moving  over  the  elliptic  surface  versus 
the  flat  plate.  The  integration  of  the  variable  convective 
phase  factor  over  the  airfoil  and  the  wake  leads  to  significant 
unsteady  effects  in  the  actual  kernel  function.  Typical  plots 
of  the  kernel  for  a  15%  thick  elliptic  airfoil  with  k  =  1 
and  Re  =  500  are  shown  in  Figs.  3.3a,  b  and  c  for  a  load 
point  at  the  leading  edge,  quarter  chord  and  mid- chord.  The 
variations  of  the  kernel  near  the  origin  due  to  finite 
Reynolds  number  do  not  appear  in  a  plot  based  on  the  scale  of 
the  airfoil  chord.  The  viscous  singular  behavior  is  illustrated 
in  the  steady  state  kernel  in  Fig.  2.3.  The  dashed  line  in 
Fig.  3.3  is  the  kernel  function  for  the  flat  plate.  The 
difference  between  the  flat  plate  and  thickness  kernel  is 
significant,  in  particular  near  the  trailing  edge  and  in  the 
wake.  While  these  kernel  variations  are  primarily  due  to 
finite  thickness  the  actual  steady  and  unsteady  load  distribu¬ 
tions  presented  in  the  following  section  depend  on  interactive 
viscous  and  thickness  effects  in  a  complicated  way  that  is 
not  fully  understood.  The  competing  effects  of  viscosity  and 
geometry  are  known  experimentally  (see  Sections  C  and  D,  below) 
and  it  is  one  of  the  important  triumphs  of  viscous  thin 
airfoil  theory  that  these  effects  can  be  calculated. 


C.  Steady  State  Results 

The  kernel  function  for  the  steady  state  problem  is 
identical  in  form  to  the  flat  plate  viscous  kernel;  i.e., 
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F1gure  3.1b.  -  S^oftt^urface  speeb  on  a 


Joukowski 
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al  convective  phas 
liptic  airfoil  (x 


Figure  3.3b.  -  Load  point  at  quarter  chord  (y 


Figure  3.3c.  -  Load  point  at  mid-chord  (y  =  0). 


(see  Ref.  2) 


K(X)  =  k  Unlxl  +  e**  K0(S|*|)3  (3-23) 


where  x  is  the  normalized  potential  coordinate  and  a  is  the 
Reynolds  number  based  on  the  radius  (approximately  quarter 
chord)  of  the  airfoil  circle  in  the  w-plane  (see  discussion  in 
Section  III.B).  The  load  distribution  on  an  airfoil  with 
finite  thickness  can  be  expressed  in  terms  of  the  corresponding 
flat  plate  solution  in  appropriate  coordinates.  From  (3.6) 
we  have 


-  %  ■  Sfr)  (3.24) 

KCO  Vco 

where  J^x)  satisfies  the  flat  plate  integral  equation 
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f  .2%)  K(x-y)dy  -  -  5T(x)  (3.25) 
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and 
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9s 


(3.26) 


is  the  modified  upwash  with  thickness  effects  included.  For 
a  thin  airfoil  at  constant  angle  of  attack,  a  ,  we  have 


3^*) 


a 

1  +  (h’)2 


(3.27) 


where  h’  is  the  local  slope  of  the  thickness  distribution 
and  must  be  expressed  in  terms  of  the  potential  coordinate  x  . 

From  our  previous  flat  plate  results  we  know  that  the 
solution  of  (3.25)  has  a  square  root  singularity  at  the  edges, 
although  the  trailing  edge  singularity  is  significantly 
weaker  than  the  leading  edge  singularity  for  moderately  large 
Reynolds  numbers.  It  follows  from  (3.24)  that  the  actual 
load  distribution  tends  to  zero  for  any  edge  that  is. 
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not  mathematically  sharp  (i.e.,  that  does  not  have  a  cusp). 

This  point  is  illustrated  in  Fig.  3.4a,  b  and  c,  where  the 
lift  distribution  is  given  for  a  Joukowski  airfoil  with 
cusped  trailing  edge  flying  forward  and  backward,  and  for 
a  near  elliptic  cross  section.  Note  the  tendency  of  the 
load  distribution  to  become  singular  at  the  trailing  edge  in 
Fig.  3,4a.  Because  of  the  cusped  edge  the  effect  of  viscosity 
is  exaggerated  and  the  overall  lift  and  moment  are  over¬ 
estimated. 

The  summary  plots  of  the  lift  curve  slope  in  Figs.  3.5 
and  3.6  illustrate  the  competing  effects  of  viscosity  and 
geometry  more  clearly.  Several  important  points  concerning 
thickness  are  evident  from  Fig.  3.5.  First  we  note  the 
direct  comparison  of  viscous  and  potential  theory  results  for 
the  standard  Joukowski  airfoil  with  cusped  trailing  edge. 

Except  for  differences  of  0(1/ /Re)  the  results  are  the  same. 
The  increase  in  lift  curve  slope  with  increasing  thickness 
predicted  by  potential  theory  with  the  Kutta  condition  is 
not  realized  in  practice  for  actual  airfoils.  The  experimental 
data  shown  in  Fig.  3.5  are  for  the  NACA  4  and  5  digit  series 
airfoils  (Ref.  17).  The  lift  curve  slope  is  less  than  the 
flat  plate  airfoil  value  and  even  decreases  with  increasing 
thickness  ratio.  For  the  NACA  63  through  66  series  (Ref.  17, 
p.  130)  there  is  a  slight  increase  with  thickness  but  the 
actual  values  are  typically  closer  to  the  value  predicted  by 
flat  plate  thin  airfoil  theory  (2tr  per  radian  or  0.1097  per 
degree) .  The  competing  effects  of  viscosity  and  geometric 
thickness  lead  to  a  nearly  universal  value  for  the  lift  curve 
slope  on  thin  airfoils,  that  operate  at  high  Reynolds  number. 
This  is  probably  the  reason  that  flat  plate  airfoil  theory 
has  served  the  unsteady  aerodynamicist  so  well. 

The  final  curve  in  Fig.  3.5  is  the  lift  curve  slope  for 
an  elliptic  airfoil  section.  The  elliptic  cross  section  is 
considered  to  be  a  reasonable  model  of  an  actual  airfoil  with 
the  displacement  thickness  added  on.  The  trend  with  thickness 
is  very  similar  to  the  trend  of  the  4  and  5  digit  series 
data.  The  actual  levels  are  higher  than  should  be  calculated 
by  a  factor  of  0(1/ /Re)  (3  or  470  in  the  present  calcula¬ 
tions).  Refinements  could  be  introduced  by  an  improved 
treatment  of  mean  boundary  layer  effects.  It  must  be 
emphasized,  however,  that  the  correct  trend,  due  to  geometric 
thickness,  is  calculated  with  only  the  inclusion  of  the 
"direct  viscous  effect"  on  a  modified  body  shape  (the  elliptic 
section  that  emulates  the  thickness  boundary  layer).  The 
improvement  over  the  results  of  potential  theory  with  the 
Kutta  condition  is  significant. 

Our  final  set  of  steady  state  results  is  shown  in  Fig.  3.6. 
The  lift  curve  slope  is  plotted  versus  the  modified  eccentricity 
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Fiqure  3.5.  -  Lift  curve  slope  versus  thickness  ratio.  Comparison  of  theory 
and  experimental  results  for  the  NACA  4  digit  series  airfoils 
(Ref.  17,  p.  130). 


parameter  e'  that  is  linearly  related  to  e  in  the  Joukowski 
mapping  [see  Eq.  (3,17)]  (i.e.,  e'  =  “e/^max)*  For  e*  “  "If 
we  obtain  the  conventional  Joukowski  airfoil  with  cusped  trail¬ 
ing  edge.  For  e'  =  +1  the  leading  edge  has  a  cusp  and  for 
e'  =  0  the  section  is  elliptic.  The  trend  of  the  calculated 
results  is  in  approximate  agreement  with  experimental  data  of 
Smith  (Ref.  28)  shown  in  Fig.  3.7.  In  particular,  we  call 
attention  to  the  difference  in  lift  curve  slope  for  the  cusped 
trailing  edge  and  the  elliptic  section.  The  calculated  re¬ 
duction  in  CL  is  15%  compared  to  a  measured  207,  reduction. 

Again  the  overall  theoretical  levels  are  higher  than  measured 
values  for  the  cusped  edge  because  of  exaggerated  viscous 
effects.  The  cusped  leading  edge  airfoil  has  a  further  re¬ 
duction  in  lift  curve  slope  of  about  10%  from  the  elliptic 
section.  The  experimental  data  of  Smith  for  a  NACA  0012 
section  flying  backwards  is  shown  for  comparison  in  Fig.  3.7. 
The  cusped  leading  edge  has  a  strong  non-linear  behavior  for 
small  angle  of  attack  and  it  is  difficult  to  infer  a  zero 
angle  of  attack  slope.  Our  best  estimate,  indicated  in 
Fig.  3.7  is  about  50%  of  the  value  for  the  same  airfoil 
flying  forward.  The  theoretically  calculated  reduction  in 
Fig.  3.6  is  about  25%.  It  is  interesting  that  for  larger 
angles  of  attack  the  NACA  0012  flying  forward  and  backward 
has  a  lift  curve  slope  that  is  almost  the  same.  The  non¬ 
linear  leading  edge  aerodynamics  are  beyond  the  scope  of  the 
linearized  viscous  thin  airfoil  theory. 


D.  Oscillating  Airfoils  -  Comparison  with  Experiment 

The  calculations  presented  in  this  section  are  obtained 
by  direct  numerical  solution  of  the  complete  integral  equation 
of  viscous  thin  airfoil  theory  [see  Eq.  (3.3)]  with  the 
kernel  function  defined  by  Eq.  (3.5).  The  basic  procedure 
is  discussed  in  Section  III. A.  For  the  purpose  of  comparison 
with  recent  experimental  data  presented  in  Ref.  24,  we 
consider  the  case  of  ri^id  body  pitch.  The  unsteady  downwash 
function  [see  Eq.  ( 3 . 6 ) j  for  the  pitching  airfoil  is: 


3T- 


1  i  k(x  -  Xp) 

.1  +  (h')2  Q0 [1  +  (h1 ) 2 ]^  _ 


(3.28) 


where  x  ranges  from  -1,  at  the  leading  edge,  to  + 1  at  the 
trailing  edge.  It  is  the  normalized  chordwise  distance  in_ 
the  physical  airfoil  plane  and  like  the  local  slope  h1  ,  it  ^ 
is  a  function  of  the  potential  coordinate  x  through  the  mapping 
defined  by  Eqs .  (3.16)  and  (3.17).  The  parameter  Xo  rs 
location  of  the  pitch  axis.  For  pitch  about  quarter  chord 
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X0  =  -0.5  .  Calculations  are  presented  for  the  elliptic  air¬ 
foil  and  the  Joukowski  airfoil  with  cusped  trailing  edge 
(t  =  127o) .  Twenty-five  Chebyshev  polynomials  are  used  to  repre¬ 
sent  the  load  distribution  and  the  computational  Reynolds  number 
is  500.  Extensive  calculations  have  shown  that  the  numerical 
results  are  sensibly  converged  with  errors  of  order  (1//Re). 

A  typical  steady  state  convergence  plot  for  the  flat  plate 
airfoil  is  shown  in  Fig.  2.4,  p.  18. 

Typical  chordwise  load  distributions  for  the  elliptic 
airfoil  and  the  Joukowski  airfoil  with  cusped  trailing  edge 
are  shown  in  Fig.  3.8a  and  b  for  a  reduced  frequency  of  0.2. 

The  pitch  axis  is  at  the  quarter  chord.  These  results  may 
be  compared  qualitatively  with  the  upper  surface  pressure 
measurements  of  Davis  and  Malcolm  (Ref.  24)  presented  in 
Fig.  3.9.  A  direct  quantitative  comparison  cannot  be  made 
because  the  Mach  number  is  0.5  in  the  experiment  versus 
Mach  zero  in  all  of  our  calculations.  However,  there  is  a 
strong  similarity  between  the  measured  and  calculated  results. 
The  important  point  to  note  is  the  reduced  loading  on  the 
elliptic  cross  section  over  the  entire  chord,  in  particular 
the  out  of  phase  load  over  the  aft  portion  of  the  airfoil. 

The  out  of  phase  components  of  the  overall  lift  and  moment  are 
quite  sensitive  to  geometric  shape,  as  we  shall  see  below. 

We  reiterate  a  point  made  previously  that  the  elliptic  cross 
section  is  considered  to  be  more  representative  of  the  NACA 
64A010  airfoil  with  boundary  layer  displacement  thickness 
added  on  than  the  Joukowski  airfoil  with  cusped  trailing 
edge.  The  cusped  edge  is  directly  responsible  for  the 
exaggerated  aft  section  loading  in  viscous  thin  airfoil 
theory.  A  direct  application  of  the  Kutta  condition  to 
potential  theory  leads  to  a  similar  exaggeration  of  the 
loads.  Even  though  the  experimental  pressure  (see  Fig.  3.9) 
over  the  aft  20%  of  the  airfoil  chord  is  almost  completely 
recovered,  the  application  of  the  "least  singular"  condition 
on  the  potential  flow  pressure  distribution  at  the  trailing 
edge  will  lead  to  an  overall  estimation  of  the  loads  that  is 
much  too  large.  This  point  was  most  clearly  brought  out  in 
the  case  of  steady  potential  flow  with  Kutta  condition 
for  airfoils  with  finite  thickness  (see  Fig.  3.5). 

The  most  impressive  display  of  unsteady  calculations  is 
presented  in  Fig.  3.10a,  b,  c  and  d,  and  should  be  compared 
with  the  corresponding  experimental  results  presented  in 
Fig.  3.11a,  b,  c,  and  d,  from  Ref.  24.  The  airfoil  is 
oscillating  about  the  quarter  chord  at  reduced  frequencies 
between  zero  and  0.3.  The  in-phase  and  out-of-phase  components 
of  the  lift  and  the  moment  about  the  leading  edge  are  plotted 
as  a  function  of  reduced  frequency.  Results  are  presented 
for  the  Joukowski  airfoil  with  cusped  trailing  edge,  the 
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Figure  3.9.  -  Upper-surface  unsteady  pres 
subsonic  theory.  (Ref.  10) 
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Figure  3.11.  -  Experimental  results  for  64A010  airfoil.  (Ref.  24) 
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elliptic  cross  section  and  the  flat  plate  airfoil  with 
Kutta  condition.  The  experimental  results  are  for  the  NACA 
64A010  airfoil  at  M  =  0.5  and  Reynolds  number  between 
2.5  and  10  million.  The  theoretical  curves  shown  in  Fig.  3.11 
were  calculated  in  Ref.  24  by  solving  Possio's  integral 
equation. 

The  overall  qualitative  features  of  incompressible 
viscous  thin  airfoil  theory,  inviscid  thin  airfoil  theory 
with  Kutta  condition  and  the  experimental  results  at  M  =  0.5 
are  the  same.  Important  quantitative  differences  are  noted 
upon  closer  examination.  For  example,  the  in-phase  component 
of  the  lift  (Fig.  3.10a)  on  the  elliptic  section  is  about 
15°/o  less  than  the  lift  on  the  Joukowski  airfoil.  A  similar 
reduction  in  the  magnitude  of  the  leading  edge  moment  is 
noted  in  Fig.  3.10c.  Note  also  that  for  reduced  frequencies 
greater  than  0.05  the  lift  and  moment  on  the  ellipse  are  less 
than  the  corresponding  values  for  the  flat  plate  airfoil  with 
Kutta  condition.  If  we  examine  the  experimental  results 
presented  in  Fig.  3.11a  and  c  closely,  we  note  first  of  all 
that  the  symbols  denoting  experimental  results  are  approxi¬ 
mately  157,  of  the  magnitude  of  the  data  plotted  and  further¬ 
more  that  in  each  graph  the  theoretical  curve  tends  to  the 
high  side  of  the  lift  data  and  the  low  side  of  the  moment 
data.  The  reduction  in  the  lift  and  moment  that  we  have 
calculated  for  the  elliptic  cross  section  are  of  the  same 
order  of  magnitude. 

The  out-of-phase  components  of  the  lift  and  moment  are 
plotted  in  Fig.  3.10b  and  d  and  should  be  compared  with  the 
corresponding  experimental  and  theoretical  results  in 
Fig.  3.11b  and  d.  In  general  the  out-of-phase  loading  is 
much  more  sensitive  to  geometry  and  reduced  frequency  than 
the  corresponding  in-phase  loading.  The  viscous  thin  airfoil 
results  for  the  cusped  trailing  edge  are  remarkably  close  to 
the  flat  plate  results.  Such  detailed  comparison  is 
considered  to  be  somewhat  fortuitous.  However,  we  note  that 
the  elliptic  section,  that  has  a  smaller  positive  loading 
over  the  aft  section  (see  Fig.  3.8a),  has  an  overall  lift 
coefficient  that  is  more  negative  and  a  smaller  moment  about 
the  leading  edge.  Also,  in  units  of  reduced  frequency,  the 
zero  crossing  of  the  lift  and  moment  for  the  elliptic  section 
is  approximately  0.05  greater  than  the  "sharp"  trailing  edge 
results.  A  comparison  of  theory  and  experiment  in  Fig.  3.11b 
and  d  reveals  that  the  zero  crossing  of  the  experimental 
results  is  also  displaced  to  the  right  of  flat  plate  (Possio) 
theory  by  about  0.05  in  units  of  k  . 

The  improvement  obtained  in  the  out-of-phase  component 
of  the  lift  and  moment  that  we  have  calculated  for  the  ellipse 
is  considered  to  be  the  most  significant  result,  to  date,  of 
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viscous  thin  airfoil  theory.  We  concur  with  the  discussion 
by  the  authors  of  Ref.  24  that  these  global  unsteady  effects 
are  in  large  measure  due  to  non-uniform  convection  of 
disturbances.  To  the  extent  that  the  loading  is  primarily 
incompressible  in  origin  the  important  phase  lag  is  in  the 
convection  time  for  a  vortex  shed  from  a  load  point  on  the 
body  to  some  point  downstream.  (Note  the  discussion  of  the 
phase  factor  in  Section  II. B,  p.  14.)  The  vortex  is  in 
general  accelerated  over  the  body  and  retarded  in  the  near 
wake.  If  compressibility  is  included  in  the  theory,  there 
will  be  additional  phase  lags  due  to  differences  in  upstream 
and  downstream  propagation  speeds. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  theory  of  incompressible  two-dimensional  viscous 
thin  airfoils  has  been  formulated  in  detail  and  numerical 
calculations  have  been  presented  for  the  symmetric  Joukowski 
family  of  airfoil  sections  in  steady  flow  and  undergoing 
simple  harmonic  motion.  Comparisons  have  been  made  with 
experimental  data  for  steady  and  unsteady  flow.  The  main 
conclusions  of  this  work  are  summarized  below: 

1.  Viscous  thin  airfoil  theory  can  be  reduced 
to  the  solution  of  an  integral  equation  that 
is  formally  the  same  as  the  familiar  equation 
of  inviscid  flat  plate  airfoil  theory.  The 
effect  of  viscosity  is  to  reduce  the  order  of 
the  singularity  in  the  kernel  function  and  to 
establish  a  unique  solution  for  airfoils  of 
arbitrary  cross  section  without  an  auxilliary 
uniqueness  criterion  such  as  the  Kutta 
condition. 

2.  Steady  flow  viscous  thin  airfoil  calculations 
for  a  cusped  trailing  edge  airfoil  with  finite 
thickness  are  in  agreement  with  the  results  of 
potential  flow  theory  with  Kutta  condition. 

The  calculated  slope  of  the  lift  curve  is  much 
too  large  and  has  the  wrong  trend  with 
increasing  thickness  ratio  when  compared  with 
experimental  data.  Any  practical  attempt  to 
include  geometric  thickness  effects  in  steady 

or  unsteady  theory  must  account  for  the  effective 
increased  thickness  of  the  aft  section  of  the 
airfoil  due  to  mean  boundary  layer  displacement. 
With  potential  theory,  the  Kutta  condition  must 
at  least  be  modified  and  with  viscous  thin  air¬ 
foil  theory  an  effective  trailing  edge  bluntness 
must  be  accounted  for. 
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3.  The  "direct"  effect  of  viscosity  ’n  the 
Stokes  layer  over  the  entire  airi  ’.1  surface 
is  the  physical  origin  of  circulatory  lift 
for  any  Reynolds  number. 

4.  The  calculated  out-of-phase  loading  on  an 
oscillating  elliptic  airfoil  is  substantially 
different  from  the  corresponding  loading  on 

a  Joukowski  airfoil  with  cusped  trailing  edge 
or  the  flat  plate  airfoil.  The  differences 
are  due  to  a  substantial  reduction  in  the  aft 
section  loading  and  agree  quantitatively  with 
the  differences  between  experimental  data 
and  direct  solutions  of  Possio's  integral 
equation  at  a  Mach  number  of  0.5. 

It  is  recommended  that  viscous  thin  airfoil  theory  be 
extended  in  two  directions  of  practical  importance. 

1.  Incorporate  viscous  and  geometric  effects 
simultaneously  in  three-dimensional  wing 
theory . 

2.  Incorporate  viscous  and  geometric  effects 

in  Possio's  integral  equation  for  compressible 
flow  and  extend  the  theory  towards  the  tran¬ 
sonic  regime. 
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with  kernel  function  whose  solution  is  obtained  with  a  modal  expansion  tech¬ 
nique  familiar  from  flat  plate  thin  airfoil  theory.  The  effect  of  viscosity 
is  to  change  the  order  of  the  singularity  in  the  kernel  function  such  that  a 
unique  solution  is  obtained  for  any  cross  sectional  geometry  without  using  an 
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auxilliary  uniqueness  criteria  like  the  Kutta  condition  or  principle  of 
minimum  singularity^  The  principal  unsteady  effect  related  to  thickness  is 
"''fo" introduce  an  explicit  phase  factor  in  the  kernel  function  that  is 
proportional  to  the  time  for  vortex  disturbances  to  be  convected  from  the 
load  point  to  the  upwash  point. 

"  T'Viscous  thin  airfoil  calculations  for  an  airfoil  with  sharp  trailing  edge 
are  in  agreement  with  the  results  of  potential  theory  with  Kutta  condition. 

The  overall  steady  loads  calculated  with  either  theory  are  in  general  much  \ 
larger  than  measured  values  and  have  the  wrong  trend  with  increasing  ' 

thickness  ratio.  The  competing  effects  of  thickness  (both  geometric  and 
boundary  layer  displacement)  and  viscosity  for  an  actual  airfoil  are  such 
that  the  measured  loads  are  in  better  agreement  with  flat  plate  theory  than 
with  the  theoretical  results  for  finite  thickness  airfoils  with  sharp  trailing 
edges,..., _ - — - — - - - -  “  - - - ' 

iscous  thin  airfoil  steady  and  unsteady  calculations  for  an  airfoil 
with  elliptic  cross  section  are  in  much  better  agreement  with  experimental 
results^  The  slope  of  the  steady  state  lift  curve  decreases  with  increasing 
>tfnckness  ratio  in  agreement  with  the  experimental  trend.  The  out  of  phase 
'loads  calculated  for  an  oscillating  elliptic  section  are  significantly 
different  from  either  flat  plate  results  or  those  for  the  Joukowski  airfoil 
with  sharp  trailing  edge.  The  magnitude  of  the  difference  is  in  quantitative 
agreement  with  the  difference  between  calculations  based  on  Poisso's  integral 
equation  and  recent  experimental  data. 

It  is  concluded  that  viscous  thin  airfoil  theory  is  a  practical  tool 
for  introducing  simultaneously  the  effects  of  viscosity  and  geometric 
thickness  in  two-dimensional' unsteady  aerodynamic  theory. 
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